Quantum tunneling of two coupled single-molecular magnets 
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Two single-molecule magnets are coupled antiferromagnetically to form a supramolecule dimer. 
We study the coupling effect and tunneling process by means of the numerical exact diagonalization 
method, and apply them to the recently synthesized supramoleculer dimer [Mn^. The model 
parameters are calculated for the dimer based on the tunneling process. The absence of tunneling 
at zero field and sweeping rate effect on the step height in the hysterisis loops are understood very 
well in this theory. 

PACS numbers: 75.45+j, 75.50Xx. 
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Nanometer-sized magnetic particles and clusters have 
generated continuous interests as study of their proper- 
ties has proved to. ba^ientifically and technologically 
very challenging .H'H'EluHfl Up to now magnetic molecu- 
lar clusters have been the most promising candidates for 
observing quantum phenomena since they have a well de- 
fined structure rjith well characterized spin ground state 
and anisotropy.Q One of the well studied systems is single- 
molecule magneUXSMM) [MaiOaClCChCCHa^dbm^] 
(short for Mn 4 ).ffll9El The molecule has well isolated 
ground state with a half integer spin S = 9/2, and 
magnetization tunneling is observed at zero magnetic 
field. Very recently a supramolecular dimer of two .SMMs 
[Mn4] 2 was reported to be synthesized successfully.ll3 The 
antiferromagnetic coupling between two SMMs leads to 
this dimer with a spin singlet ground state and makes 
the quantum tunneling quite different from SMMs M114. 
The coupling also makes this dimer an excellent candi- 
date for studying quantum tunneling in a system of two 
truly and coupled identical particles. Quantum tunnel- 
ing of magnetization can be advantage for some applica- 
tions of SMMs in providing quantum superpositions of 
states required for quantum computing.EjE-3 So the cou- 
pling effect in quantum tunneling between SMMs is a 
very important issue in application of integrated molec- 
ular magnets. In this paper we first study the tunneling 
process in one SMM Mn4 with spin 5 = 9/2. A local 
stray field has to be introduced to-explain the tunneling 
of the ground state at zero ficld.B Then we study the 
coupling effect of two SMMs and observe a novel triangle 
tunneling process. We apply our observation to the newly 
synthesized supramolecular dimer of two SMM [Mn4] 2 
and make use of the two triangle tunneling processes to 
deduce the the model parameters from the experimental 
data, and, furthermore, explain the sweeping rate effect 
in the derivatives of hysteresis loops. 

We first start with a biaxial model for a SMM Mn4 
with spin 5 = 9/2, 



Hi = -DSi 



(B + h) (1) 



B = Be z is the external magnetic field along the z 
axis. The term g/igfioSi ■ h is the local stray field in- 
teraction between the SMMs and the environment. For 
simplicity we denote the energy eigenstate of the biax- 
ial model I to) by its dominant S z — m component, and 
to = —S, —S + 1, • • • ,S. The E term and the stray field 
may lead to some minor correction to these states. If the 
stray field is not included, it is well known that for a half 
integer spin the tunneling between the st ate |n -j >) and \S) 
is quenched due to the parity symmetry.! 1 5 tt L6 H 1 1 1t can be 
proved simply that, for any integer n, we always have 
(S\ (Hi) n \M) = with M = -S + 1, -S + 3, • • • , S. 
The equality indicates that there is no connection or no 
tunneling occurs between these states \—S) and \M) . Ex- 
perimentally quantum tunneling was observed in a SMM 
Mn4 between the states \—S) and \S)B So there must 
be a small internal magnetic field by the nuclear spins of 
the Mn ions and/or the dipole-dipole interaction between 
different molecules. We model the interaction as a local 
stray field h with a random Gaussian distribution with 
the equal width a in tbaee directions as we did for the 
molecular magnets Fegjia 



P(h) 



1 



(27TC7 2 ) 3 / 2 



exp[-h 2 /2cr 2 ]. 



(2) 



where i = 1 or 2 referring to the two SMMs in the 
dimer, D and E are the axial anisotropic constants, 



A transverse component of such a field may lead n to a 
tunneling splitting at zero field as observed in Rem. In 
this paper we take tho-parameters for a SMM Mn4 D = 
0.762K, E = 0.0317KJ3E3 and a = 0.035tJB The result- 
ing tunneling splitting for the ground states | — 9/2) and 
|9/2) at zero field is V( A o) = 3.280x 1(T 7 K, and that for 
the states |-9/2) and 1 7/2) is y/{Af) = 1.52627 x 1(T 5 K 
by using the exact diagonalization method jHiherc (• • • ) 
stands for the averaging over the stray field.Ej Thus the 
local stray field may cause a tunneling splitting between 
the ground states. 

Following Wernsdorfer et al., the two SMM Mn4S in the 
dimer [Mn4] 2 are coupled via a weak antiferromagnetic 
superexchange coupling J. Thus the model Hamiltonian 
for the dimer is 



H = ^ + H 2 + JSi ■ S 2 



(3) 
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FIG. 1: Spin state energy of the model Hamiltonian (Eq. 
(H)) for [Mn4]2 without local stray fields as a function of ap- 
plied magnetic field by taking D = 0.762K, E = 0.0317K, 
and J = 0.1K. The two triangles are related to the tun- 
neling (1) from |-9/2, -9/2) + to |-9/2,9/2) + at point 1, 
then to |9/2, 9/2) + at point 3, and (2) from |— 9/2, -9/2) + 
to |-9/2, 7/2) + at point 2, then to 1 7/2, 7/2) + at point 4. 
The Point 5 is from |-9/2,9/2) + to |9/2,7/2) + . The reso- 
nance fields for the five points are -0.335T, 0.233T, 0.335T, 
0.861T, and 0.943T, respectively. 



where Si = S2 — 9/2. For each dimer there are (2Si + 
1)(2S 2 + 1) = 100 energy eigenstates. Like in a SMM 
Mn<4, each state can be labelled approximately by two 
predominant quantum numbers |mi,m 2 ) for two SMMs 
with mi. 2 = -9/2, -7/2, • • • , 9/2. Without the coupling 
J the states \mi,m2) and \m%,mx) are degenerated. As 
the two SMMs can be regarded as truly identical particles 
there exists the permutation symmetry between particle 
1 and 2 and the eigenstates may possess parity symme- 
try. Thus the eigenstates for the system are denoted by 
l m ii TO 2) + for even parity and |mi,m2)_ for odd parity. 
The antiferromagnetic coupling J may remove the de- 
generacy of these two states, but the parity in the states 
remain unchanged. Even when we take into account the 
the coupling J and the transverse terms the states be- 
come the linear combination of all possible states, for 
simplicity, we still use the two dominant quantum num- 
bers to represent the states. All the energy eigenvalues 
by neglecting the local stray fields are plotted in Fig. 1. 
The average over the local stray field does not move the 
position of energy level crossing. 

Before explaining the experimental observation from 
the dimer we first consider the effect caused by the cou- 
pling J between the two particles. Assume the tunneling 
between the states |mi) and \my) under a sweeping field 
B( = — c;,t) and the tunneling splitting between the two 
states is A, the pair of the splitting energy eigenvalues 
near the resonant point can be written as 




2 

with Cb = g/iBl^afidB /dt, and the two states are given by 
4> l ±{t) = {±c±\m l )+c T \m v ))/V2 (5) 

with c± = yli (mi — rni>)ci,t/ \J~ (mi — mi>)cbt\ 2 + A 2 . 

Before the resonant tunneling, the initial state is at 
I mi) , i.e., at t = —00, 4>+(t) — * \mi') an d <f>-(t) ~ * 
|mi) ; after the tunneling, at t = +00, cf> + (t) — ► |mi) 
and 4>-{t) — > —\my). When two identical parti- 
cles are put together there are four possible states: 
|+, +) + = 4>+ ® 4> 2 + with the energy 2e + , |+, — ) ± 
= {4>\ <8> 4>^ ± </>! ® </>+) fy/2 with the energy e + + e_, 
and | — , — ), = (j>l_ ® cf) 2 _ with the energy 2e_. We denote 
the even and odd parity of the states by the subscripts 
±, respectively. The energies vary with time t, and 
are plotted in Fig. 2, for an illustration, by choosing 
mi = m-2 = —9/2 and mi' = my — 7/2 for the model 
in Eq.(^). The tunneling splitting between the two 
states |+, +) , and | — , — ), is 2 A, the double of a single 
particle as expected. To see the coupling effect of two 
particles, we plot the energy eigenvalues for several 
different couplings in Fig. 2. It is obtained by the exact 
diagonalization of 100 x 100 matrix for the Hamiltonian. 
The two states |+, — ) ± are degenerated for J = 0.0. A 
very little coupling J may remove the degeneracy of the 
two states. The state | + , — )_ has odd parity and does 
not take part in the tunneling process as other three 
states have even parity. It is shown that the coupling J 
leads to two consequences: (1) The tunneling splitting 
from |— ,— ), and |+,+), decreases very quickly with 
increasing J, and almost closes for J > 0.3 x 10~ 5 K. In 
the dimer of [Mn^ the coupling J « 0.1K and tends to 
suppress the tunneling at this point completely. (2) The 
tunneling splitting from | — , — ) + to |+, +) + occurs at two 
separated points via a intermediate state |+, — ) , . The 
coupling J provides an inner bias field to expel the two 
resonant points away from the original ones of |+,+), 
and |+, — ) , . This triangle process reflects the structure 
of the tunneling of the two identical particles. In the 
language of mi and 7712, the process from |mi,mi), to 
|mi/,mi')i is described as follows: tie first resonant 
tunneling occurs from |mi,mi), to |mi,mi/) + , and 
the second one follows from |mi,mi') , to |mi/,mi/) + . 
The explicit tunneling from |mi,mi) + to |mi',mi'), is 
suppressed completely by the coupling J. 

Now we are ready to analyze the quantum tunneling 
in the dimer. Typical hysteresis loops in magnetization 
versus sweeping external field applied along the easy axis 
are observed. These loops display step-like features sepa- 
rated by plateaus. The step heights become temperature- 
independent below 0.3K, but depend on the sweeping 
rate of magnetic field c = dB/dt. Derivatives of the loops 
at different sweeping rate reflect that quantum tunneling 
occurs at several points, but is absent at zero field. At 
high field the initial state is |— 9/2,— 9/2) which has 
even parity. Due to the permutation symmetry of iden- 
tical particles all the tunneling to the states |mi,m.2)_ 
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FIG. 2: Energy levels of the states |— 9/2, -9/2> + , |7/2, 7/2) + 
and |— 9/2, 7/2) ± versus sweeping field at different coupling 
J of two coupled identical particles as an example to illus- 
trate the tunnenling process from |mi,rai) + to |m 1 /,m 1 /) + 
via |mi,m 1 /) + . (x-axis: Energy/10 _4 K, y-axis: Magnetic 
Field/5xl(T 6 T) 



with odd parity in this system are prohibited. The tun- 
neling process in the dimer can be understood essentially 
by two triangle processes as shown in Fig. 1. Starting 
from the initial state | — 9/2, —9/2),, the first level cross- 
ing happens at magnetic field b\ = — 0.33K, which is 
from |-9/2, -9/2) + to |-9/2,9/2) + at point 1, and the 
dual resonant point is from |— 9/2, 9/2) + to |9/2,9/2) + 
at point 3 63 = +0.33K in the first triangle process. The 
energy of the intermediate state | — 9/2, 9/2) + is indepen- 
dent of the external field. The resonant field for the 
points 1 and 3 are 61,3 = ±9J/2g[iBfi<o from the model 
(Eq.(3)). Thus J is calculated to be O.IK as Wernsdorfer 
et al. did. The finite coupling does not lead to a tunnel- 
ing splitting at this point, which can be proved explicitly: 
for an integer n we have 

+ (-9/2, -9/2| H n |-9/2, 9/2) + = (6) 

when the stray field is absent. The tunneling splitting at 
points 1 and 3 are caused by the local stray field. These 
two points were consistent with Ref.E3. The second pro- 
cess is from | — 9/2, — 9/2) + to 1 7/2, 7/2) + via an interme- 
diate state I — 9/2, 7/2) , . The two energy level crossings 
are from |— 9/2, — 9/2) , to |— 9/2, 7/2> + at point 2, and 
from |-9/2, 7/2) + to [7/2, 7/2) + at point 4. We take 
the parameters for D and E for a SMM M114, and find 
that the calculated resonant fields are 62 = 0.233T and 
64 = 0.861T-™ which is very closed to the experimental 
data 0.87T.E3 As the points 2 and 3 are very close such 
that the resonant peaks are smeared to a broaden one. In 
Wernsdorfer et al's paper they neglect the transition from 
1—9/2, 7/2> + to |7/2,7/2) .and thought that the third 
peak in the Fig. 4 of Ref.E3 are caused by those from 
|-9/2, -9/2) + to 1-9/2, 5/2) + and from |— 9/2, 9/2) + to 
|7/2,9/2) + (i.e. the point 5 in Fig. 1). D was calculated 
to be 0.72T by neglecting the transverse component E. 
After a detailed analysis we found that the transverse 
component E remove the degeneracy of | — 9/2 , 9/2) , and 



|-9/2, 9/2)_ . Even if we take D = 0.72T we find that 
b 2 = 0.20T, 64 = 0.833T and b 5 = 0.9138T which is 
larger than 0.87T. If we fix the point 5 at 0.87T, D is cal- 
culated to be 0.664T smaller than 0.72T. The tunneling 
from |-9/2, -9/2) + to |-9/2,5/2) + belongs to another 
the triangle process and the splitting which is caused by 
the stray field is much smaller than those at points 2 
and 4. On the other hand we anticipate that the weak 
coupling between two SMMs does not affect the intrinsic 
properties of a SSM in the dimer too much. Our calcula- 
tion shows that it is reasonable to take the parameters D 
and E from the SMM Mn4 for the model of the dimer of 
the two SMMs with a weak coupling. It is worth point- 
ing out that the coupling J can also drive the tunneling 
splitting between some states such as | — 9/2 , +7/2) + and 
|9/2,7/2) + . However these tunnelines. do not contribute 
significantly to what observed in RefJia We do not discuss 
them here. 

After determining the positions of the resonant points 
and model parameters we are in a position to calculate 
the tunneling splitting, which determines the transition 
rate in the Landau-Zener model. The exact diagonaliza- 
tion method is applied to calculate the energy eigenvalues 
at different external field. The sampling average is taken 
for the local stray field. For each sampling we calculate 
the energy levels as in Fig. 1 and find the energy split- 
ting A„ at each resonant point. More than 1000 sampling 
arc taken to calculate the averaging tunneling splitting 
yj (A^ ) for each distribution width a. The calculated tun- 
neling splittings are listed in Table I. 

Table I: The calculated tunneling splitting yj (A 2 ) in 
unit 10 _5 K at different resonant points with different dis- 
tribution width a of the stray field h by using the exact 
diagonalization method. ( D = 0.762K, E = 0.0317K.) 
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2.81552 


1.52671 


0.010 


0.01960 


2.19227 


0.01960 


2.83557 


1.59807 


0.020 


0.03207 


2.19487 


0.03207 


2.84268 


1.59723 


0.035 


0.04687 


2.20155 


0.04687 


2.87249 


1.61907 


0.050 


0.06264 


2.78338 


0.06264 


2.98306 


1.68816 



The derivatives of the hysteresis, loops at different 
sweeping rates in Fig. 4 of Ref.tHl indicate that the 
peak heights in the derivatives depend on the sweep- 
ing rate. The height of the first peak decreases with 
the increasing rate while oppositely the second peak in- 
creases. This phenomenon can be understood qualita- 
tively in the modified Landau-Zener model. In princi- 
ple the time-evolution of the spin system can be reached 
by solving the 100(= (2S + 1)(2S + 1)) coupled time- 
dependent Schrodinger equation. As the tunneling split- 
ting is very small, the coupled equations near the two 
resonant states can be reduced to an effective two-level 
system with the Hamiltonian 

H = ( (m v +m 2 ')(c b t + h z ) A(h)/2 \ 

e// I A(h)/2 ( mi +m 2 )(c b t + h z ) ) 
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FIG. 3: Calculated derivative of the hysterisis loop at different 
sweeping fields based on the modified Landau-Zener model. 



rate is proportional to the reverse of the sweeping rate 
c = dB/dt, approximately. The larger the sweep- 
ing rate is, the less particles tunnel to the new state. 
The step height is related to the transition rate by 
AM = PLzJ2i( m i ~ m 'i)- The presence of the local 
stray field will smear the "iamp" of the magnetization 
around the resonant point £3 At a field b around the 
resonant point, the variation of the magnetization due 
to quantum tunneling can be approximately given by 
M(b) ~ AM dh 2 P(h 2 ),which leads to the deriva- 
tive of the hysteresis loop around the resonant point, 
dM/db ~ AMP (b). The calculated results are plotted 
in Fig. 3. Comparing with Fig. 4 in Ref.L3, we find that 
the numerical results based a Landau-Zener model are 
consistent with the experimental observation, and essen- 
tially reflect the sweeping rate effect on the peak height 
of derivatives of hysteresis loops. 



The tunneling splitting A(h) between two states 
l m ij m 2) + an d {mi/, rri2>) + can be obtained by diagonal- 
izing the Hamiltonian with a specific field h. The state 
evolves with time, 



®eff(t) = exp 



H eff (t)dt 



®eff(t = -00), 



and the magnetization varying with time is given by 
M(t) = ($ e/ f \ Sf + S| |$ e / /) ■ The average over the stray 
field h is taken for (dM{t)/dt) = J dhP(h)dM(t)/dt. 
Physically, with the local stray field, the Landau-Zener 
transition formula is given by 

Plz = 1 - (exp[-7rA^ m ,/z/ mm /]) w it ( A^ m , ) /v mm > , 



where v mm > = 2g/i s /i 



Z)i=i2( m »' - m i) dB/dt. The 



In conclusion, we study the coupling effect of two truly 
identical particles, and analyze the quantum tunneling of 
magnetization in the supramolecular dimer of two Mn4S. 
The exchange coupling between two SMMs provides a 
biased field to expel the tunneling to two new resonant 
points via an intermediate state, and direct tunneling is 
prohibited. Based on the analysis we deduce the model 
parameters from the experimental data, and find out that 
the coupling does not change the model parameters for a 
SMM too much. Finally we point out that the sweeping 
rate effect in the derivatives of hysteresis loops can be 
explained quantitatively in the modified Landau-Zener 
model. 
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